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Abstract 


This  report  consists  of  three  parts.  In  Part  1, 
the  eilconal  and  continuity  equations  are  derived  for 
a sound  wave  propagating  through  a fluid  in  which 
there  is  shear  motion  of  low  Mach  number  and  scale 
large  compared  to  the  sound  wave  length.  On  the  basis 
of  a comparison  between  the  eikonal  equation  and  the 
Hamilton- Jacobi  equation  for  a charged  particle  in  a 
magnetic  field  it  is  shown  that  the  acoustic  ray  paths 
are  identical  with  the  trajectories  of  the  particles 
if  the  magnetic  field  is  proportional  to  the  vorticity. 

In  Part  2,  the  attenuation  and  fluctuation  in 
intensity  of  a sound  beam  scattered  in  a turbulent  med- 
ium are  expressed  in  terms  of  correlation  products  of 
the  turbulent  flow.  Expressions  are  also  derived  for 
the  phase  and  intensity  fluctuations  associated  with  the 
propagation  of  sound  through  large  scale  turbulence,  as 
discussed  in  Part  1. 

In  Part  3,  the  differential  cross  section  for 
scattering  of  sound  by  turbulence  is  developed  in  a form 
suited  to  treatment  of  anisotropic  turbulence.  A simple 
anisotropic  distribution  involving  symmetry  about  a pre- 
ferred vorticity  axis  Is  discussed. 
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Notation 


Both  vector  and  tensor  notations  are  employed 
according  to  convenience.  Vectors  are  indicated  by 
a wavy  underline.  Repeated  tensor  indices  are  to  be 
summed.  Sub-  or  superscripts  in  parentheses  are  not 
tensor  indices.  The  anticommuting  symbol 
has  the  value  one  for  <*- , is  antisymmetric 
to  permutation  of  any  pair  of  indices,  and  vanishes  if 
two  indices  are  equal.  In  terms  of  this  symbol  the 
cross  product  of  the  vectors  and  is  Cr  ~ 

Symbols  of  the  form  Indicate  integration  over 
three  dimensional  space.  The  symbol  indicates  in- 
tegration over  all  angles. 

Space  derivatives  such  as  3f/3x  * , 3 f *0,d'x/r  are 
frequently  denoted  by  and  respectively. 


On  the  Propagation  of  Sound 
in  a Turbulent  Fluid 


Part  1.  Propagation  of  Sound  through  Large  Scale  Turbulence 

1.1  When  a sound  wave  propagates  through  a fluid  in  shear 

motion,  energy  is  scattered  out  of  the  beam^-  if  the  spatial 
power  spectrum  of  the  shear  motion,  or  turbulence,  has  appreci- 
able strength  at  wave  numbers  which  are  the  same  order  of  size 
as  that  of  the  sound  wave.  When  the  dominant  eddy  size  in  the 
shear  flow  is  considerably  larger  than  the  sound  wave  length  such 
scattering  becomes  very  small,  but  there  remains  a coherent  type 
of  perturbation  which  results  in  gradual  distortion  of  the  wave 
fronts  without  loss  of  energy.  This  distortion  can  be  described 
by  a WKB  type  approximation  procedure  which  yields  an  eikonal 
equation  determining  the  phase  fronts  and  a continuity  equation 
describing  the  amplitude  fluctuations  associated  with  the  phase 
distortion. 


In  the  presence  of  a shear  velocity  field  of  low  Mach 
number  the  wave  equation  for  sound  propagation  has  the  form^ 
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(1.1) 


where  V is  the  excess  pressure,  pe  the  mean  pressure,  M the 
longitudinal  or  sound  velocity  obeying 


u. 


= 0 


(1.2) 


and  w"  the  shear  or  transverse  velocity  obeying 


t & 
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(1.3) 


The  velocity  of  sound  in  the  undisturbed  medium  is  c. 
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In  the  absence  of  the  shear  velocity  v/,  the  pressure 
and  velocity  for  a weak  plane  sound  wave  of  arbitrary  form  may 
be  taken  as 


-p- Af  C*-n«x*)  U*  = d A n,  -f  (i-hv-x*)  (l.l^.) 


where  A is  a dimensionless  constant  which  we  shall  call  the  amp- 
litude. We  now  try  to  satisfy  (1.1),  under  the  assumption  that 
the  Mach  number  of  the  shear  flow  is  small  and  that  & varies  but 
little  in  a sound  wavelength,  by  replacing  A with  the  space 
function  A(z>  and  with  the  function  ■&(£),  The  usefulness 

of  this  approach  is  that  because  varies  slowly,  the  variations 
in  A and  from  their  unperturbed  forms  will  also  be  slow. 

By  making  use  of  the  divergence  and  curl  conditions 
(1.3)  and  (1.2),  the  right  side  of  (1.1)  may  be  written 
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( V oi^ \a/c< J/3  ) 


Since,  because  of  the  low  Mach  number  associated  with  w,  this  term 
is  a small  perturbation,  we  may  replace  by  the  unperturbed 

value  within  the  limits  of  the  approximations  to 

which  (1.1)  is  valid. 

We  shall  now  substitute  into  (1.1)  the  functions 


-p  =p6  A(x)f  U^-d  AC*)#,*  f ( i-tfC*)) 


where  ■$,<*  stands  for  2-&/dx*  . The  result  Is 


(1.5) 


Af  vMf  (i,6) 


where  the  prime  denotes  differentiation  with  respect  to  the 
argument.  The  quantities  are  small  quantities 

in  this  equation  according  to  the  assumption  of  slow  variation 
of  the  shear  velocity.  Hence,  neglecting  terms  bilinear  in  these 
quantities  we  obtain  the  approximate  equation 
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Af  Af ' "+  V' f 'A  + 2 A,* „ f ' = 2 <r'[w< Af {U7) 


Vie  have  also  neglected  the  term  in  V2A  , which  will  be  small 
compared  to  A,<*  * Since  f is  a quite  arbitrary  function,  the 
coefficients  of  f"  and  f ' must  be  independently  zero;  we 
therefore  obtain  the  two  equations 


zA^S^+AV1#  =-2*"i2,<i9//rw*//fA 


(1.8) 

(1.9) 


The  first  of  these  relations  is  an  eikonal  equation  determining 
the  phase  as  a function  of  space.  Rewriting  the  second  in  the 
form 


A1),  + yzi!i  + z<t~'-&JaL$JfWuJ/a  — O 


(1.10) 


it  is  seen  to  be  a continuity  equation  relating  square  of  amp- 
litude to  phase.  Eq.  (1.10)  shows  that  the  variations  in  amp- 
litude arise  from  two  causes:  focusing  effects  associated  with 

the  curvature  of  the  wavefronts;  and  the  necessity  for  mainten- 
ance of  constant  energy  flow  through  regions  in  which  the  local 
velocity  of  propagation  varies  in  the  direction  of  propagation. 

1.2  In  this  section  it  will  be  shown  that  the  effect  on  the 

propagation  of  sound  rays  of  large  eddy  3ize  turbulence,  or,  in 
general,  shear  motion  which  varies  but  little  In  the  space  at  a 
sound  wave  length,  bears  a direct  analogy  to  the  effect  of  mag- 
netic fields  on  the  paths  of  charged  particles.  Since  the  latter 
has  been  extensively  studied  in  connection  with  fields  such  as 
cosmic  ray  research  and  particle  accelerator  design,  the  analogy 
possibly  could  be  a valuable  tool  in  visualizing  sound  propagation 
phenomena.  The  analogy  is  based  on  a similarity  between  the 
eikonal  equation  for  sound  waves  In  a slowly  moving  medium  and 
the  Hamilton-Jacobi  equation  for  the  motion  of  charged  particles 
In  a magnetic  field. 

Consider  a particle  whose  motion  is  determined  by  the 
Hamilton  equations 

^ H 

^ = r«=  d.ii) 


where  q«  is  the  position  vector  of  the  particle,  the  time  rate 
of  change  of  this  vector.  The  Hamiltonian  H and  momenta  70*  are 
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related  to  a time-independent  Lagrangian  function  L Cq*jqa)  by 


3 L/Z%  H Cr>«j<}<)=  £ V*  9k  - L (1.12) 


The  possible  motions  of  such  a particle  corresponding 
to  a total  energy  E may  then  be  found,  according  to  the  Hamilton 
Jacobi  procedure’,  by  solving  for  S the  equation 

H ( ^ £ (1.13) 


in  which  25/59*  replaces  p>*  wherever  it  appears  in  the  Hamiltonian. 

The  paths  along  which  the  particles  of  energy  E travel 
may  be  described  very  simply  in  terms  of  the  function  S.  Accord- 
ing to  Hamilton's  principle  of  least  action^  the  particle  paths 
determined  by  (1.13)  are  those  for  which  the  action  integral 

t 

(1.14) 

is  extremal  (maximum  or  minimum)  between  two  arbitrary  fixed  end 
points  q'"  and  q^.  Regarding  the  time  as  a parameter  along  the 
path,  we  may  write  q* <J£  = <1 9*  so  that,  substituting  55/39*  for 
the  extremal  condition  on  (l.lij.)  may  be  expressed  as 

y(Z) 

If  now  a wave  pattern  were  constructed  such  that  S were  every- 
where proportional  to  the  total  phase  of  the  waves,  (1.15)  would 
determine  those  paths  between  the  two  fixed  end  points  along 
which  the  total  change  of  phase  would  be  extremal.  It  is  therefore 
a statement  of  Fermat's  principle,  and  the  paths  so  defined  are 
simply  the  ray  paths  of  the  wave  system  as  ordinarily  understood. 

We  have  therefore  obtained  the  general  result  that  the 
ray  paths  of  a wave  system  are  identical  with  a family  of  trajec*- 
tories  of  particles  obeying  a Hamilton- Jacobi  equation  identical 
with  the  eikonal  equation  for  the  waves.  This  fact,  which  is 
the  basis  of  the  correspondence  between  classical  and  quantum 
mechanics,  takes  a simple  and  well  known  form  in  the  case  where 
the  Hamiltonian  contains  no  velocity  dependent  forces.  Then, 


(3)  See,  e.g.,  Goldstein,  ’’Classical  Mechanics"  Addison-Wesley, 
Cambridge  (195>1) 


(Ij.)  See  reference  (3)  above 
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the  particle  velocity  cj*  is  proportional  to  the  canonical 
momentum  73*  so  that  the  ray  paths  are  normal  to  the  phase 
fronts,  S = constant,  and  the  ray  direction  at  any  point  is 
the  direction  of  the  vector 

The  Lagrangian  of  a non-relativistic  particle  of 
mass  m and  charge  e in  an  electromagnetic  field  with  vector 
potential  A *67)  and  scalar  potential  VC$_)  is 

4.  A.1  - (1.16) 

2 <i0 

where  c0  is  the  velocity  of  light. 

The  corresponding  momenta  and  Hamiltonian  are 

Hipt^)  = C A«)V  Z + «V 


Accordingly  the  Hamilton- Jacobi  equation  is 
/ 3 S e * \ 1 , 

( Jj.  ~ /l"  * ev  = £ 

where  S is  the  total  energy  of  the  particle. 

The  approximate  eikonal  equation  for  a sound  wave 
propagating  in  a medium  undergoing  shear  flow  was  found  to  be 

9-,*  + 2.  c-'xP,*  w*  ~ / = 0 (1#8) 


(1.17) 


(1.18) 


according  to  section  1.1  Here  w$/<!  is  the  total  phase,  « is  the 
frequency,  t the  velocity  of  sound  in  the  still  medium,  and  w* 
the  3hear  velocity  field  obeying 


Since  (1.3)  is  accurate  only  to  terms  linear  in  w*  , it  may  be 
written  to  equally  good  approximation  as 


(j&j  * -+  c ~ ' W*  ) 1 — 


/ 


(1.19) 
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which,  it  may  be  verified,  differs  from  (1.3)  by  only  the  term 
cr*|*lx. 

In  this  form,  the  similarity  to  the  Hamilton- Jacobi 
equation  (1.18)  with  V-0  is  readily  apparent.  The  similarity 
is  made  more  complete  by  the  fact  that  the  vector  potential 
obeys  A when  V=constant,  in  analogy  to  (1.3)  • The  exact 
identification  between  (l.l8)and  (1.19)  may  be  shown  by  dividing 
(1.18)  by  E,  and  noting  that  -2. is  simply  the  square  of  the 
ordinary  momentum  of  the  particle,  which  we  may  denote  by  P*. 

In  terms  of  this  quantity  (1.18)  may  be  written 

(S,*  - eq0~'A«  )*/  P 1 = I (1.20) 


which  is  identical  with  (1.19)  if 

IP  » i£j0<  - e CT'A*  IP  = (1.21) 


Thus,  paths  of  sound  rays  in  a fluid  in  which  the  local  velocity 
is  W<  are  identical  with  the  trajectories  of  particles  of  charge 
e and  momentum  P moving  in  an  electromagnetic  field  corresponding 
to  the  vector  potential  A*  = da  6~‘Pw*/g, 

Several  limitations  on  the  analogy  must  be  kept  In  mind. 
First  we  have  assumed  that  the  velocity  field  w*  changes  but  little 
in  a distance  of  a sound  wavelength,  and  is  small  compared  to  c. 
Secondly  we  have  not  as  yet  Included  the  effects  of  time  variation 
of ; our  discussion  of  the  Hamilton-Jacobi  equation  was  based 
on  the  assumption  of  time-independent  potentials.  Actually,  the 
requirements  of  large  eddy  size  and  small  Mach  number  imply  that 
the  time-variations  in  w*  will  be  very  slow  indeed  compared  to 
those  of  the  sound  wave  and  except  for  special  cases  they  will 
not  result  in  significant  departure  from  the  time  independent 
equations. 


The  magnetic  field  corresponding  to  the  vector  potential 
in  (1.21)  is 


as  Cur-I  A.  = — ( do  C*'P/0  CO I 
eH/c.P  = - t-'Vf  (1.22) 


where  X is  the  vorticity  associated  wlthw.  Hence  the  results 
of  this  section  may  be  summed  up  by  saying  that,  within  the 
limitations  noted  above,  the  ray  paths  of  a sound  wave  in  a 
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moving  fluid  follow  the  trajectories  of  charged  particles  in 
a magnetic  field  everywhere  proportional  to  and  in  the  opposite 
direction  from  the  vorticity  vector,  the  proportionality  con- 
stant being  determined  in  terms  of  the  (arbitrary)  particle 
momentum  by  (1.22). 5 

It  was  noted  above  that  every  turbulent  velocity  field 
w corresponds  to  a possible  magnetic  field  in  free  space,  and 
the  converse  is  also  kinetically  possible.  Hence,  the  propaga- 
tion of  sound  in  a moving  atmosphere  should  exhibit  all  the 
focussing  effects,  shadow  effects,  and  other  phenomena  encountered 
with  particles  travelling  in  magnetic  fields.  It  would  be  of 
interest  to  examine  to  what  extent  the  analogy  holds  true  when 
eddy  sizes  the  order  of  a wavelength  or  smaller  are  involved. 

This  would  correspond,  if  the  analogy  held,  to  the  quantum 
mechanics  of  charged  particles  in  magnetic  fields.  Also  of  in- 
terest is  the  question  of  whether  such  time -dependent  phenomena 
as  the  betatron  acceleration  mechanism  have  counterparts  in  the 
sound  propagation  case.  To  investigate  this  it  is  necessary  to 
extend  the  treatment  above  to  include  time  dependent  potentials. 
There  could  possibly  result  an  analogy  between  the  ray  paths  for 
the  time  dependent  case  and  the  trajectories  of  the  4“dimenai°nal 
Hamilton- Jacobi  equation  for  a relativistic  particle. 


Part  Attenuation  and  Intensity  Fluctuations 

2.1  When  a sound  wave  propagates  through  a turbulent  fluid 

medium,  the  local  intensity  of  the  wave  is  altered  both  by 
scattering  of  energy  out  of  the  beam,  as  discussed  by  the  author 
in  a previous  paper, 1 and  by  distortion  of  the  shape  of  the  wave- 
fronts  by  large  scale  shear  flow  as  indicated  in  the  previous 
section.  In  this  section,  the  intensity  variations  to  be  ex- 
pected from  both  effects  will  be  expressed  in  terms  of  the 
statistical  correlation  products  associated  with  the  turbulence. 

The  calculations  will  be  confined  to  cases  in  which  the  character- 
istic time  periods  of  the  turbulence  are  very  large  compared  to 
the  period  of  the  sound  wave;  this  represents  a realistic  re- 
striction for  the  discussion  of  sound  propagation  in  the  atmosphere 
or  sea,  where  very  small  Mach  numbers  are  ordinarily  encountered. 

In  reference  (1)  expressions  were  derived  for  the 
differential  scattering  cro3s-section  for  sound  waves  traversing 
a region  in  turbulent  motion.  It  follows  from  the  analysis 
developed  there  that  when  the  Mach  number  of  the  turbulence  is 


(5)  Note  that  there  is  no  direct  relation  between  the  speed 

of  propagation  of  the  rays  and  the  speed  of  the  particles . 
Only  the  spatial  forms  of  the  paths  are  compared-. 
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small,  the  frequency  spectrum  of  the  differential  scattering 
cross-section  is  given  by 


I 


— C<6  (0.  'Q  ) 


|UJ 


“"*)l 


(2.1) 


Here  c is  the  velocity  of  sound,  h.  a unit  vector  in  the  direc- 
tion of  the  incident  sound  wave,  »il>a  unit  vector  in  a direction 
of  scattering,  w.  the  frequency  of  the  incident  sound-wave,  and 
Wc&oj)  the  !(.- dimensional  Fourier  transform  of  the  turbulent 
velocity  field.  ICri"  &'*> ) is  related  to  the  differential  scatter- 
ing cross-section  by 


. <t)  ( I / . (h  (ii 


(2.2) 


where  T is  the  period  of  time  over  which  the  time  Fourier  analysis 
is  performed. 

The  principal  effects  on  the  scattering  of  the  time 
variation  of  the  turbulence  as  reflected  in  the  Jp-dimensional 
transform  WC&«)  are  first,  a slight  broadening  in  frequency  of 
the  scattered  wave,  and  second,  an  effect  on  the  amount  of 
scattering  at  very  small  angles  in  radians  of  the  order  of  mag- 
nitude of  the  ratio  of  the  characteristic  frequencies  of  the 
turbulence  to  w..  These  effects  will  have  a negligible  influence 
on  the  attenuation  of  a sound  beam,  which  depends  on  the  total 
cro3S-section  and  hence  we  shall  express  the  scattering  cross- 
section  in  the  time  independent  form 


<rcd':dv)  = 


n 

c1 


C2J  t , ll> 

g • W'cfc.ta  - 


a ]) 


(2.3) 


where  and  V/C-U  the  3-dimensional  Fourier  transform  of 

the  shear  velocity  is  defined  by 


W t&) 


-if  - *‘A-x  , 

Un)  aJw  c«;e  ~ d 


(2.1p) 


Equation  (2.3)  can  be  derived  by  performing  the  analysis  of 
reference  (1)  under  the  initial  assumption  that  the  shear  veloc- 
ity is  independent  of  time.  In  using  this  expression  to  calculate 
the  attenuation  of  a sound  wave  over  long  distances,  it  is  assumed, 
of  course,  that  the  turbulence  is  statistically  stationary  over 
the  regions  and  times  involved. 
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The  scattering  cross-section  as  it  appears  in  (2,3) 
may  be  expressed  readily  in  terms  of  the  correlation  products 
usually  employed  in  describing  turbulence.  Rewriting  the  ex- 
pression in  WCk), 


,(|)  JIK*  .<*■>  I.'1’* 


•w  . . 


If  the  correlation  tensor  of  the  turbulence  is  denoted  by 


4/3  C ) ~ ''v*-  ^ ~ ) f >:  + A .) 


(2.6) 


where  the  overbar  indicates  averaging  over  a suitable  region 
of  space,  the  Fourier  transform  of 


) - Cztcj  2/^CA  )e 


(2.7) 


differs  from  ) only  in  that  it  is  normalized  by  the 

total  volume  of  turbulence  considered.  Hence,  the  scattering 
cross-section  per  average  cubic  centimeter  of  turbulence  will 
be 


- . <n 


r(n 


d "> 


21  u*  / 

C *3-  ■ 


‘*\2  <*’  -T- 


«/* 


CM, 


r 0 7 \ 

lH  - £ ]) 


(2.8) 


The  total  scattering  cross -section  is  found  by  in- 
tegrating (2.8)  over  all  scattering  directions  rP': 


0/ 

> 


7t  , 4 ^ /.  | (5i  2.  (27  -r-  . (Z)  (7)  v | 2-  ,^2.7 

~ r*  **  ^ rv  ~ (2.9) 


If  the  initial  intensity  of  the  beam  of  sound  is  Ie,  the  mean 
intensity  after  traversing  a distance  / will  be 


I 


— ^ (tf 

,-rca  yy 


(2.10) 


Formula  (2.10)  will  describe  the  attenuation  correctly 
only  if  y is  greater  than  the  distance  over  which  there  exists 
appreciable  correlation  among  the  eddy  sizes  (wave-numbers ) 
contributing  to  the  scattering.  Within  shorter  distances,  there 
will  be  coherence  and  interference  between  the  waves  scattered 
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from  different  points.  The  restriction  to  distances  longer 
than  the  correlation  length  is  implicit  in  (2,9)  because  the 
space  averaging  involved  in  the  correlation  tensor  Is 

customarily  carried  out,  in  the  case  of  extended  turbulence, 
over  regions  large  compared  with  the  values  of  \ of  interest. 

The  Mach  number  of  the  turbulence  encountered  in  the 
atmosphere,  the  extended  medium  of  principal  interest,  is 
sufficiently  low  that  the  amount  of  attenuation  in  a distance 
the  order  of  several  correlation  lengths  may  be  expected  to  be 
quite  small.  In  this  case,  the  relative  fluctuations  to  be 
expected  in  the  attenuation  will  be  very  nearly  averaged  out 
in  a distance  of  transmission  long  enough  to  result  in  sig- 
nificant weakening  of  the  sound  beam.  An  estimate  of  the 
intensity  fluctuations  may  be  obtained  as  follows.  Let  the 
linear  distance  over  which  the  space  averaging  in  (2.6)  has 
been  carried  out  be  L.  Then  the  mean  scattering  cross-section 
of  a region  of  the  turbulence  one  centimeter  square  and  ex- 
tending: a distance  L in  the  direction  of  propagation  will  be 
Lf  Ch111).  The  fluctuations  from  the  mean  in  the  scattering  from 
this  "region  will  then  be  given  by  LA  a-  cnl,>)  where 


3 *.n< 


. O I U) 


(2.11) 


The  corresponding  mean-square  fluctuation  is 


<U<rc ntt"/ 

y{^J:/*CkeL ) AT^x C *. L hu,/- r>"  1 (2.12) 


where  the  averaging  is  over  a representative  distribution  of 
regions  of  linear  size  L. 


If  the  fluctuations  in  <r  are  taken  into  account,  the 
intensity  after  traversing  a distance  y will  be 


X 


Jdy 


(2.13) 


Hence,  the  mean-square  fluctuation  in  the  logarithmic  decrement 
of  the  intensity  will  be 

( /e9  Y*  )*>  = *y)1} 


(2.llf) 
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It  has  been  assumed  that  correlation  in  the  turbulence  is 
negligible  outside  the  distance  L.  Hence  if  y is  considerably 
larger  than  L we  may  expect  that  the  total  fluctuation  appearing 
in  (2.1I4.)  is  the  result  of  random  addition  of  fluctuations  from 
y/L  independent  regions.  According  to  a well-known  law  we 
should  then  have 


< O §°)*)  -({)<[*  e-CG!1') 


(2.15) 


so  that  the  ratio  of  the  fluctuation  in  logarithmic  decrement 
to  the  mean  decrement  is 


iK,  (A  (l0/i))*->  /(/*3  /<^c  (2.l6) 


It  is  apparent  from  (2.12)  that  the  calculation  of  the 
fluctuation  < [a <rcc";j j*)  requires  a knowledge  of  the  fourth  order 
correlation  product 

< * '«/*  C M ) -h  (2.17) 


which  may  be  expressed  in  terms  of  velocity  correlation  products 
by  the  use  of  (2.7)  and  (2.6).  The  evaluation  of  such  fourth 
order  correlation  functions  is  difficult  with  existing  theoretical 
and  empirical  knowledge  of  turbulence. 

2.2  As  shown  in  section  1.1  the  phase  of  a sound  wave  propa- 

gating In  a fluid  in  which  there  is  low  Mach  number  turbulence, 
the  scale  of  which  is  large  compared  to  the  sound  wavelength, 
becomes  distorted  according  to  the  elkonal  equation 

^*^+24''^  W*  - / =0  (1  8) 


where,  if  w Is  the  frequency  of  the  sound  wave,  is  the  total 

phase.  The  shear  velocity  field  is  w^,  If  the  unperturbed 
sound  is  a plane  wave  propagating  in  the  direction  of  the  unit 
vector  the  unperturbed  phase  will  be  - £■*  so  that, 
writing  -&■  — $-+i9-i0)  , the  eikonal  equation  may  be  written 


+ 2 c~'  0«w/*  + w*  ) = 0 (2.18) 


If  the  distortion  Is  gradual,  corresponding  to  the  assumption  of 
low  Mach  number,  <f>  will  vary  quite  slowly  as  a function  of  space 
and  hence  the  terms  $,«.  and  will  be  small  compared  to 
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the  other  terms  in  (2.18).  Therefore,  the  phase  distortion 
will  be  determined  to  a first  approximation  by  the  equation 


+ <"r\«  Wo < = o 


(2.19) 


If  the  value  of  the  phase  is  given  on  some  wavefront  by,  say, 
^=o,  then  according  to  (2.19)  the  phase  further  along  the 
direction  of  propagation  will  be  given  by 


4>  C * + 2-n)  = - c"' 


,x'+ni 

. 'Wen;  c*;  d.  cx-n  > 


(2.20) 


where  is  a point  on  the  initial  wavefront,  is  the 

component  of  w in  the  direction  of  n.  and  the  line  integral  is 
over  a length  L along  the  direction  of  n . Although  (2.19)  is 
approximately  correct  everywhere  in  the  region  of  shear  motion, 
it  must  be  noted  that  the  integrated  form  (2.20)  will  be  correct 
only  when  the  change  of  direction  of  the  wavefronts  indicated 
by^cx)  is  small;  i.e.,  the  error  introduced  in  passing  from 
(2.18)  to  (2.19)  is  essentially  that  the  initial  direction  of 
propagation  £,  appears  in  (2.20)  instead  ot  the  local  direction, 
which  is  the  resultant  of  the  initial  direction  and  the  total 
distortion  which  has  occurred  up  to  the  point  of  interest. 

According  to  (2.20)  the  mean  square  phase  deviation  at 
the  point  x - x°  + £! i.  is 


<1  <^0,  (*.')  w(n,  Cx  " ) > <J  <£'-*)  d C*"-h)  (2,21) 


where  the  averaging  is  over  a statistically  representative 
ensemble  of  shear  flow  systems  so  as  to  provide  a valid  measure 
of  the  phase  distortion  in  the  turbulent  flow.  The  quantity 
<wc Cjs/ivv/o,,  Cx".>  is  just  a diagonal  component  of  the  two  point 
correlation  tensor  <w .ex')  If  the  statistics  of  the 

turbulence  are  spatially  homogeneous,0  this  tensor  is  a function 
of  x'-x"  only  and  reduces  to  the  single  argument  correlation 
tensor  defined  by  (2.6).  In  this  case  (2.21)  reduces  to 

<L4c*>ix>  n*  Cx'-x")  (*'■£.)  <JCx"-£0  (2.22) 


(6)  This  does  not  imply  isotropy 
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This  expression  is  easily  reduced  to  single  integrals.  If 
we  imagine  a plane  of  integration  in  which  y~  2L  •£  and  h~;2A"'2 
are  perpendicular  cartesian  coordinates,  the  integration  in 
(2.22)  is  over  a square  of  side  L , oriented  with  the  y and  » 
axes.  It  is  readily  seen  by  drawing  the  diagonal  lines 
y~2=A  that  the  weight  with  which  the  value  A occurs  in  the 
square  is  L-Ui  . Since  A ranges  from  L to-L  and  since 

> sr  the  integral  in  (2.22)  may  be  re- 

written 

<£iz 5c>°+Lh)Jl>  = /eL^CAn  )(L-AjdA  (2.23) 


where  we  have  noted  in  writing.  that  the  integration 

is  along  n . 

When  the  distance  of  propagation  L is  considerably 
longer  than  the  maximum  distance  over  which  there  exists 
correlation  in  the  turbulence,  CXr\ ) will  be  negligible  for 
values  of  A greater  than  L and  hence-,  for  this  case, 7 


M)  ( 2.2ij.) 


Under  the  assumption  that  L is  large,  the  second  term  will  be 
small  relative  to  the  first  and  hence  the  plausible  result 
that  the  root-mean-square  phase  deviation  varies  asVL,  and  as 
a quantity  associated  with  the  mean  velocity  fluctuation  length 
and  magnitude. 

According  to  section  1.1  the  local  intensity  of  the 
distorted  wave  (more  accurately,  the  square  of  the  amplitude) 
is  related  to  the  phase  deviation  by  the  continuity  equation 

= -V*t£  - 2 v^,  (1  o) 


where  e is  the  logarithm  of  the  intensity.  Under  the  approxima- 
tion employed  in  deriving  the  phase  deviation  we  may  replace 
(1.9)  by  the  approximate  form 

€'«n*  - ac'1  n*  h*  (2.25) 


(7)  A similar  derivation  was  carried  out  for  a specific  model 
of  turbulent  flow  by  Blocklntzevi  "Acoustics  of  an 
Inhomogeneous  Moving  Medium."  Translation  by  R.  T,  Beyer 
and  D.  Mintzer  under  RAG  Contract  N7-0NR-35808  (August  1952) 


i4- 


so  that,  following  the  analysis  just  given. 


£ Q ) 


X 


LV'0Cj 


(2.2b) 


will  represent  the  approximate  logarithmic  intensity  increment 
after  propagating  a distance  L into  the  moving  medium.  Writing 
out  Vx  in  the  form  3 ■*,)*+  sVo  xt)1  where  n lies 

along  the  x?  axis,  and  using  expression  (2.20)  for  fie*), 


vV<*)  = 


.*/  3wjC?j 
c — — 

c>  Xj 


(2.27) 


Hence, 

< (2>28) 


The  term  Cv/jO^+in)—  vy, cx6))/ d depends  only  upon  variations,  in 
the  direction  of  propagation,  of  w,,  the  component  of  shear 
velocity  in  the  direction  of  propagation.  Its  physical  sig- 
nificance is  that  it  represents  the  change  in  amplitude  required 
to  maintain  a constant  transfer  of  energy  as  the  velocity  of 
propagation  along  a line  of  propagation  varies  according  to  w^c*; 
It  is  apparent  that  this  term  is  of  a fluctuating  nature,  and 
in  view  of  the  spatial  homogeneity  of  the  turbulence  which  we 
have  assumed,  it  will  not  make  any  cumulative  contribution  to 
€ Cxv  L£),  and  will  not  be  considered  further. 

The  first  term,  Involving  derivatives  with  respect  to 
andxZj  represents  changes  in  local  intensity  due  to  fluc- 
tuations in  the  spacing  of  ray  paths  as  a result  of  the  curvature 
of  the  phase  fronts.  In  fact,  3*0/cs»x,)l-*-3V  ^ gives  just 

the  local  curvature  of  a phase  front.  This  type  of  intensity 
fluctuation  may  be  regarded  as  the  resultant  of  a series  of 
localized  focusing  and  defocusing  actions.  Let  the  integral  in 
(2,28)  be  denoted  by  eTc<)  . Then,  if  we  employ  the  notation 

L 9 /Cd*, )x-t-  K )*■  j f Cj>t)  = ) /p  on  ) 


for  compactness, 


<[€TCx6+tn)Jx>  = J-  J-  Cxi  ^XCx^y  4x^41 

y x 


/ 

* 


(2.30) 
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*7 


(2.30) 


represents  the  mean  square  intensity  fluctuation  due  to  wave 
front  curvature.  Using  ("2. 20), 


< ^ <£ ) <*')>  - t~*  j^o  / <w*- cx ) vs^cx')  > c/ x/ 


(2.31) 


I'Jriting 


<^wjI C* ) vv/cxh)  = d1f  0£-x'0 


(2.32) 


we  note  that  this  correlation  product  can  depend  only  on  the 
difference  of  the  two  coordinate  points  in  the  case  of  homogen- 
eous turbulence;  then, 


cfx  ' 


(2.33) 


where  y = > Y £.'~ £ ° , and  the  variables  of  integration  have 

been  transformed  to  give  more  convenient  limits.  The  integral 
in  (2.33)  is-  similar  to  that  in  (2.22)  except  that  the  area  of 
integration  is  a rectangle  in  this  case  instead  of  a square. 
Using;  a geometrical  construction  similar  in  principle  to  that 
employed  in  the  evaluation  of  (2.22),  one  finds,  with  3ome 
manipulation, 

y y' 

- T^fCl^Cy-y'-XU* 

'O 

It  will  be  noted  that  this  expression  is  symmetric  in  y 
and  y'.  Denoting  it  by  c-CyyO  , 


<CeT  = f^J^o-Cy,  y'jJyJy' 


(2.35) 


As  in  the  case  of  the  phase  deviation,  this  expression  for  the 
intensity  fluctuations  in  terms  of  integrals  over  second  order 
velocity  derivative  correlation  products  can  be  simplified  when 
the  propagation  distance  L is  much  longer  than  distances  over 
which  there  exists  correlation  in  the  turbulence.  For  such  long 
distances,  following  the  argument  used  in  obtaining  the  approxi- 
mation (2.24.)> 

O-  Cy,y')  — Cy+y'>/t~fe*£)JA  ~ )y~y'l 


(2.36) 
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Approximately  correct  results  will  be  obtained  by  substituting 
this  asymptotic  form  for  G-Cy>y')  into  (2.35)*  The  reason  is 
that  since  G-Cy^yO  grows  withy  and  y ' , the  contribution  to 
(2.35)  from  values  so  small  that  (2.36) is  not  valid  will  be 
minor.  Hence  setting 

F-  f0°°f  (*h)d\  (2.37) 

we  have 

<rexfxvi.h)3*>  ~^//iL£fy+yo-)y-y'i3  JyJy  ' (2>38) 

Considering  this  as  an  integral  over  a square  of  side  L , it  may 
be  handled  as  were  the  previous  cases,  except  that  for  evalua* 
tion  of  the  term  in  y+y',  the  weighting  of  the  lines  y+y'~\  will 
be  h-lL-Al , and  A will  range  from  0 to  2L . Then, 

L L <*A-2jf*o_-A>ns=  ft-3  (2.39) 

Hence, 

<£6r  ~ f (2.4.0) 

where  F is  defined  in  terms  of  the  quadratic  correlation  of  the 
velocity  derivatives  by  (2.37) , (2.3 2),  (2.29).  Thus  the  root- 
mean-square  intensity  fluctuation  will  increase  as  Lr.  The  rapid 
increase  of  the  intensity  fluctuations  relative  to  the  phase 
deviation  is  due  to  the  fact  that  the  fluctuation  at  distance  is 
the  cumulative  result  of  wavefront  curvature  along  the  entire 
path  of  propagation  as  well  a3  at  the  point  L . 


Part  3»  Scattering  by  Anisotropic  Turbulence 

In  a previous  paper^-  an  explicit  expression  was  derived 
for  the  scattering  of  sound  in  a region  containing  isotropic 
turbulence  of  arbitrary  spectrum  shape.  In  this  section,  the 
effects  of  anistropy  on  the  scattering  are  discussed. 

3.1  In  most  physical  examples  of  turbulence,  there  is  a 

significant  degree  of  anisotropy  because  the  mechanism  by  which 
energy  is  supplied  to  the  turbulence  tends  to  produce  vorticity 
in  certain  preferred  directions.  In  the  process  of  transfer  of 
the  energy  from  large  eddies  to  small  the  anistropy  becomes 
obscured  gradually,  and  near  the  upper  limit  of  the  turbulence 
spectrum  there  is  usually  a high  degree  of  isotropy.  In  terms 
of  spatial  distribution,  it  may  be  expected  that  the  maximum 
anistropy  exist  near  the  vorticity  producing  boundaries  while 
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us  the  vorticity  diffuses  away  from  the  boundary,  the  anisotropy 
become  obscured.  This  is  true,  for  example,  in  the  case  of  at- 
mospheric turbulence  resulting  from  vorticity  produced  at  the 
ground. 


In  speaking  of  the  isotropy  or  anisotropy  of  the  turbulence 
in  a region  of  fluid  we  refer  not  to  the  instantaneous  character 
of  the  flow  but  to  statistical  properties  which  are  defined  as 
averages  over  either  a long  period  of  time,  or,  as  is  practically 
equivalent  in  this  case,  a very  large  volume  of  fluid  containing 
many  such  regions  as  we  are  examining.  Alternatively,  the  statis- 
tical characteristics  may  be  expressed  as  averages  over  an  ensemble 
of  flow  system^  the  ensemble  being  so  chosen  as  to  represent  the 
fluctuations  encountered  in  time  averaging  or  large  scale  space 
averaging.  In  what  follows  the  ensemble  average  of  a quantity  A 
will  be  denoted  by  <A> . 

Consider  a shear  motion  for  which  the  velocity  is  given 
by  . The  velocity  field  may  be  described  in  terms  of  a 

Fourier  representation  by 


V/,  Ck^)  — C 2r<j'ljy*/„C£Jt)e. 


(3.1) 


The  condition  that  the  motion  be  a pure  shear  flow  is 

3 woi  C 2£  ■> 

— — o 

3 y* 


(3.2) 


or 


k*W«  = o 


(3.3) 


In  the  reference  mentioned,  it  was  shown  that  when  a 
plane  sound  wave  of  frequency  uj,  propagating  in  the  direction  of 
the  unit  vector  nc,>  is  scattered  in  a region  of  a fluid  undergoing 
such  shear  flow  at  low  Mach  number,  the  power  spectrum  of  the 
energy  scattered  per  unit  solid  angle  in  the  direction  *V11  is  given 
by 


i * 


7 rw 
< < 


U)  (.hi 


Cn-n  ) 


ID  , , 


u;  n i 

( - oo x , 


(34) 


I Cj*,' is  so  normalized  that  the  differential  scattering  cro3s- 
soction  is  given  by 


(>)  <») 


cr(h  ; n ) 


= ~rf 


X (S"‘ 


j j d u» 


(3.5) 


18 


where  ~c  is  the  interval  over  which  the  time  integration  in 
(3.1)  is  carried  out. 

In  order  to  investigate  the  effects  of  anisotropy 
on  the  scattering  predicted  by  (3*4) » Xt  Is  convenient  to 
express  the  Fourier  elements  C h,*)  in  terms  of  the  vorticity. 
The  vorticity  axial  vector  is  defined  by 

( 3 • k ) 

'V^r  (-jf)  *•*)  ~t  k/3  Wli  to) 

where  the  quantity  is  antisymmetric  to  permutation  of  any 

two  indices  and  €.)23=l.  From  this  it  follows,0  as  a result  of 
(3.3),  that 

0.7) 


Substituting  this  expression  in  (3»^-)» 

(ii  a,  troj*,  <ii  01.  xi  -1 


x cd'l a%  - ^<au-j2w>Vl€«>* 

with  k ~ 11“  In1-  ^oV"!  . 

Because  of  the  antisymetry  of  in  c<  and  cr  the  term  in 

vanishes,  leaving 


T/_c<  u>  , _ O)  u/  tf  Sx>  ur />l“lnu,~c^nt'\  u1 

— <T  <>>1  •&  ) I %-«,<*.  >V  r)# 

-»  f\  \ — . 4 


X C IcWjQ.  - 


8) 


or,  introducing  the  vector  voloclty  potential  V^CXjtj  by 
W<r  ) = 1 k/t  V«  ) 

W<r  C*, t ) = 9 *)/*  * ^ 


(3.9) 


so  that 

yCcx}t)-VlVrCxjt) 


Vr  Cjk,*)=  -A*Var  (3.10) 


(8)  This  is  readily  shown  by  multiplying  (3*8)  by  and  using 

the  relation  = 


(3.8)  becomes 
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iGfX«o  = 


6*  1 


V,( 


. U)  tij 

l<j?n  - 


) w| 


or 


(3.11) 


This  form,  or  (3*8)  displays  the  angular  behavior  of 
the  scattering  more  immediately  than  the  original  expression 

(3.11) »  and  shows  that  £tl>and  ntU  play  symmetrical  roles  in  the 
scattering  formula.  The  quantity  in  the  straight  brackets  has 
a very  simple  geometrical  significance.  It  is  equal  to  the 
volume  of  the  parallelopided  formed  by  the  vectors  G%  ti."  and 
V.  Thus  the  angular  dependence,  aside  from  the  detailed  struc- 
ture of  V,  is  the  combined  result  of  two  factors;  first,  the 
quantity  jV"* n'9 * 11  which  favors  scattering  in  forward  and  backward 
directions  while  indicating  zero  90°  scattering,  and  second,  the 
parallelopided  volume  which  is  maximum  when  the  angles  between 
the  three  vectors  are  right  angles  and  hence  indicates  zero 
scattering  at  zero?  and  180°,  and  also  in  the  case  when  either 
n‘"  or  Allies  in  the  direction  of  the  vorticity  component  for  the 
particular  wave  number  involved.  Thus  in  the  case  of  a shear 
flow  in  which  there  is  a preferred  direction  of  vorticity  it  is 
to  be  expected  that  the  total  scattering  will  be  maximum  when 
the  direction  of  the  incident  wave  is  normal  to  the  vorticity 
direction  and  that  the  scattered  radiation  will  be  most  intense 
in  or  near  the  plane  normal  to  the  vorticity.  This  specific 
angular  behavior  would  be  superimposed  on  the  angular  dependence 
associated  with  isotropic  turbulence. 

Both  the  vorticity V and  the  velocity  potential  V have 
vanishing  divergence  if  the  velocity  field  itself  is  to  be 
divergenceless.  However,  the  divergenceless  character  of  g as 
defined  in  (3*9)  *s  independent  of  any  condition  onV(,  Hence, 

(3.11)  will  describe  the  scattering  from  some  possible  incom- 
pressible flow,  irrespective  of  whether  or  not  y is  divergence- 
less. In  the  latter  case,  however,  the  vorticity  is  not  related 
to  y.  in  the  simple  fashion  indicated  in  (3*10). 

The  advantage  of  describing  the  shear  velocity  field  by 
the  inherently  divergenceless  expression  (3*9)  is  that  one  need 
not  impose  explicitly  on  the  statistical  ensembles  describing 
turbulent  flows  the  condition  (3.3)  that  the  velocity  vector  and 
wave  vector  be  perpendicular.  As  an  example,  for  which  no  ob- 
servational validity  is  claimed,  we  shall  examine  the  scattering 
from  a velocity  distribution  which  describes  a particular  type  of 
hypothetical  turbulence  with  a preferred  axis  of  vorticity  and 
which  leads  to  especially  simple  expressions.  Consider  the 


(9)  Strictly  speaking,  it  cannot  be  immediately  concluded  that 

the  scattering  at  zero  degrees  vanishes  without  ascertaining 
the  behavior  of  }/(&>*)  for  very  3mall  & . It  is  apparent, 

however,  as  was  pointed  out  in  reference  (1),  that  the 
frequency  spectrum  of  the  radiation  in  the  forward  direction 
must  be  very  sharp. 
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velocity  spectrum 

W,  (*,  w)  = ka  hyf  C & UJ  ) ■f  Ckj'^)  ( 3 .12 ) 


in  which  tf°'is  a unit  vector.  For  a fixed,  constant  this 
spectrum  describes  motion  confined  to  directions  lying  in 
the  plane  perpendicular  to  n*".  If  it  is  assumed  that  the  en-^ 
semble  distribution  of  nCM  is  independent  of  &<*>  and  that  the 
anisotropy  of  the^  turbulence  is  completely  described  by  the 
distribution  oft!*';  i.e.,  if  depends  only  on  k*j  «*> 

and  -f  C4y“>  is  uncorrelated  with  *1%  the  ensemble  average  of  the 
scattering  will  be  given  by 


WCP)5,“j)-  J ^ — =^,/w|-wb)  (3.13) 


where  F Ok)]/1*)  ~ K y 

For  thi3  situation  it  is  apparent  that  the  dependence 
of  the  scattering  on  the  direction  of  the  incident  wave,hf'i  is 
entirely  contained  in  the  factor  j*>  which  is  the 

mean  square  volume  of  the  parallelopiped  bounded  by  the  three 
unit  vectors.  This  quantity  is  easily  calculated  for  a given 
distribution  of  rf01  around  some  preferred  direction.  In  view  of 
the  simplicity  of  this  example  It  would  be  of  interest  to  see 
to  what  extent  axially  symmetric  turbulent  flow  occurring  In 
practice  could  be  fitted  to  the  distribution  described.  The 
principal  lack  of  generality  in  the  axially  symmetric  distribu- 
tions described  by  (3.12)  with  the  specializations  we  have  noted 
is  that  a definite  weighting  is  given  to  wave-vectors  * according 
to  the  angle  they  make  with  the  axial  direction  tt'»  as  may  be 
evidenced  by  forming  the  energy-spectrum  corresponding  to  (3.12): 

1 )£/ <*J  * I ~ L (-{*•  *3  I f Oj  | 1 (3, lk) 


This  weighting,  which  favors  wave -vectors  lying  in  the  preferred 
plane  of  motion,  is  not  an  unreasonable  one,  a priori,  but  it 
should  be  emphasized  that  no  empirical  validation  of  this  dis- 
tribution is  implied. 


Appendix 


In  the  analysis  of  Section  2,2  there  appeared 
several  integrals  of  the  form 


I,  — Jg  GC/x-yDdxrfy 


(A.l) 


or 


I*  - X X H 


(A.  2} 


These  expressions  are  readily  converted  to  single  Integrals 
by  first  finding  the  weight  with  which  given  values  of  the 
arguments  A'  1.x.- y I or  A=xt-y  occur  in  the  rectangular  regions 
of  integration.  Considering  Ij.  , first,  it  appears  from 
Figure  1 that  in  region  II  of  the  domain  of  integration  the 
weight  corresponding  to  the  argument  value  X is  bAA,  the 
area  of  the  parallelogram  of  base  and  height  b , Similarly, 
neglecting  areas  of  order  (dAX,  the  weights  appropriate  to 
regions  I and  III  are  (b-IAi)dA  and  Ca-AJdA  respectively.  In 
region  I,  lAI  varies  from  o to  b,  in  region  II  from  o to 
and  In  region  III  from  a.-b  to  a.  Hence  I,,  expressed  as  the 
sum  of  Integrals,  over  the  three  regions,  is  given  by 


I,  = / 6-CA)  (b-A)dA-t-  f G-(A)MA  d-  f <?  (A)  Ca-'A  )dA  (A. 3) 

0 e 


If  we  note  that 


-d-b 

J G-CMbdA  = J^CA)Ca-A)dA-J  G-CAXa-b-A;dA 


this  may  be  written  in  the  form, 


This  expression  is  symmetric  In  a.  and  b (as  is  readily  verified 
for  the  last  term)  and  hence  Is  valid  also  for  Whena^b, 

I,  = (A. 6) 


which  agrees  with  the  result  of  the  derivation  described  in  the 
text  for  the  integration  over  a square. 
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Considering  now  the  evaluation  of  X* , it  appears 
from  Figure  2 that  in  region  I the  weight  corresponding  to 
the  argument  value  is  while  in  region  II  It  is 

Cza.~X)dA  . in  region  I X ranges  from  o to  a and  in  region  II, 
from  a to  2a.  Hence,  both  these  weights  can  be  expressed  by 
the  single  function  (a.- /a- Al)<U  • Then, 

r2a 

X*  — J ^ H (A)  C<x  - la -Xi)d  A (A. 7) 


This  method  of  evaluating  Iz  can  be  extended  to  rectangular 
domains  of  integration  also. 
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